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Abstract 

The hitting time of a classical random walk (Markov chain) is the time required to detect the presence 
of - or equivalently, to find - a marked state. The hitting time of a quantum walk is subtler to define; 
in particular, it is unknown whether the detection and finding problems have the same time complexity. 
In this paper we define new Monte Carlo type classical and quantum hitting times, and we prove several 
relationships among these and the already existing Las Vegas type definitions. In particular, we show 
that for some marked state the two types of hitting time are of the same order in both the classical and 
the quantum case. 

Further, we prove that for any reversible ergodic Markov chain P, the quantum hitting time of the 
quantum analogue of P has the same order as the square root of the classical hitting time of P. We also 
investigate the (im)possibility of achieving a gap greater than quadratic using an alternative quantum 
walk. In doing so, we define a notion of reversibility for a broad class of quantum walks and show how to 
derive from any such quantum walk a classical analogue. For the special case of quantum walks built on 
reflections, we show that the hitting time of the classical analogue is exactly the square of the quantum 
walk. 

Finally, we present new quantum algorithms for the detection and finding problems. The complexities 
of both algorithms are related to the new, potentially smaller, quantum hitting times. The detection 
algorithm is based on phase estimation and is particularly simple. The finding algorithm combines a 
similar phase estimation based procedure with ideas of Tulsi from his recent theorem |Tul08| for the 2D 
grid. Extending his result, we show that for any state-transitive Markov chain with unique marked state, 
the quantum hitting time is of the same order for both the detection and finding problems. 

1 Introduction 

Many classical randomized algorithms are based on random walks, or Markov chains. Some use random 
walks to generate random samples from the Markov chain's stationary distribution, in which case the mix- 
ing time of the Markov chain is the complexity measure of interest. Others use random walks to search 
for a "marked" state in the Markov chain, in which case the hitting time is of interest. In recent years, 
researchers studying quantum walks have attempted to define natural notions of "quantum mixing time" 
[NV001 lABN+Oll I AAKVOT] and "quantum hitting time" [AKR051 I5ze04l IMNRS07] and to develop quantum 
algorithmic applications to sampling and search problems. 

Ambainis |Amb04j was the first to solve a natural problem — the "element distinctness problem" — using 
a quantum walk. Followin g this, quantum walk algorithms were discovered for triangle finding [MSS07J, 



matrix product verification BS06 , and group commutativity testing |MN07| . All of these are "hitting time" 
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applications involving quantum walk search on Johnson graphs — highly-connected graphs whose vertices 
are subsets of a fixed set and whose edges connect subsets differing by at most two elements. Quantum 
walk algorithms for the generic spatial search problem [AA05j were given by Shenvi et al. |SKW03] on the 
hypercube, and by Childs and Goldstone [CG04] and Ambainis et al. AKR05) on the torus. Szegedy |Sze04j 
showed that for any symmetric Markov chain and any subset M of marked elements, we can detect whether 
or not M is nonempty in at most (of the order of) the square-root of the classical hitting time. To achieve this 
goal, Szegedy defined the quantum analogue of any symmetric Markov chain. Later Magniez et al. [MNRS07 
extended this to define the quantum analogue of the larger class of irreducible Markov chains. 

Unresolved by Szegedy's work is the question: with what probability does the algorithm output a marked 
state, as opposed to merely detecting that M is nonempty? This issue was addressed by Magniez et 
al. [MNRS07J, who gave an algorithm which finds a marked state with constant probability but whose 
complexity may be more than the square root of the classical hitting time. Indeed, for the y/N x y/N 
grid their bound is 0(iV) whereas the classical hitting time is <d(N\ogN). The algorithms of Ambainis et 
al. }AKR05| and Szegedy jSze04j perform better on the grid if there is a unique marked state: they find the 
marked state in time 0(y/~N • log AT). (The case of multiple marked elements may be reduced to this case 
at the cost of a polylog factor in run-time.) For some time it remained unclear if one could do better, until 
Tulsi |Tul08j showed how to find a unique marked element in time 0(y/NlogN). His algorithm seems to be 
something of a departure from previous quantum walk algorithms, most of which have been analyzable as 
the product of two reflections a la the Grover algorithm [Gro96J. The 2D grid was the canonical example 
of a graph on which it was unknown how to find a marked state quantumly with the same complexity as 
detection. Tulsi's result thus raises the question: is finding ever any harder than detection? 

In this paper we address several questions related to classical and quantum hitting times. In the literature 
on Markov chains, hitting time is usually defined as the complexity of the natural Las Vegas algorithm for 
finding a marked element by running the chain. We first give an alternative definition based on the Monte 
Carlo version of the same algorithm. To our knowledge, this variant of the hitting time has not been 
considered previously. We show that for some marked state, the two hitting times are of the same order 
(Theorem E3]). 

Within the setting of abstract search algorithms presented by Ambainis et al. [AKR05] . we introduce 
quantum analogues of the two classical hitting times (Definition 13. 2p . The analogue of the Las Vegas version 
was already present in Szegedy's work |Sze04j , whereas the other is new. Unlike in the classical case, detection 
and finding are substantially different problems in quantum computing. We address both problems here. 

For the detection problem, we introduce a new algorithm Detect based on phase estimation which is 
similar to the approach of Magniez et al. MNRS07 . Our algorithm can detect the presence of a marked 
element in the starting state. The advantages of this algorithm are its simplicity and the fact that its 
complexity is related to the new Monte Carlo type quantum hitting time (Theorem 13. 5p . This is an 
improvement over the Szegedy detection algorithm whose complexity is related to the potentially larger Las 
Vegas quantum hitting time. 

We then present a variant of the above algorithm, called Rotate, which can be used for the more 
difficult problem of finding, and whose complexity is also related to the Monte Carlo type quantum hitting 
time (Theorem I3.9j) . This improves the finding algorithm due to Ambainis et al. whose complexity was 
characterized by a potentially larger quantity, the inverse of the smallest eigenphase of the search algorithm. 
Our algorithm is also simpler. Combining Rotate with the ideas in the Tulsi algorithm for the 2D grid, we 
can find a marked element with constant probability and with the same complexity as detection for a large 
class of quantum walks — the quantum analogue of state-transitive reversible ergodic Markov chains. 

As in the classical case, for some marked elements the two types of the quantum hitting time are of the 
same order (Fact l3~10l and Theorem I3.14|) . For any reversible ergodic Markov chain P, we prove that the 
quantum hitting time of the quantum analogue of P is of same order as the square root of the classical 
hitting time of P (Theorem 13. 15p . Moreover, for the Las Vegas hitting times they are exactly the same. 

Finally, we investigate the (im)possibility of achieving a greater than quadratic gap using some other 
quantum walk. For this we consider general quantum walks on the edges of an undirected graph G; these 
were defined, for example, in the survey paper of Ambainis [Amb03j . see also |San08j . We define a quite 
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natural notion of reversibility for general quantum walks. We conjecture that for any reversible quantum 
walk U2 on an undirected graph G, there exists a reversible ergodic Markov chain P on G such that for 
every marked state, the quantum hitting time of U% is at most the square root of the classical hitting time 
of P. We are able to prove this in the special case of quantum walks built on reflections f Theorem 13. 18|) . 
Our proof introduces a classical analogue of such quantum walks which might be of independent interest 
(Definition 13. 5[) . Curiously, the classical analogue is reversible if and only if the quantum walk is reversible 
(Lemma 15717)1 . 



2 Classical hitting times 

Let P be an ergodic and reversible Markov chain over state space X = {1, .. . ,n}, which we identify with 
its transition probability matrix. We suppose that the eigenvalues of P are nonnegative, by replacing P 
with (P + 1) /2 if necessary. More generally, one may also assume that the eigenvalues of P are at least a, 
where a £ [0, 1), by replacing P with ((1 — a)P + (1 + a)I)/2 if necessary. We make this further assumption 
when needed, for instance in Section 2J Let ir denote the stationary distribution of P. Let P_ z be the 
(n — 1) x (n — 1) matrix we get by deleting from P the row and column indexed by z. Similarly, for a vector 
v, we let V- Z stand for the vector obtained by omitting the z-coordinate of v. 

Claim 2.1. The eigenvalues of P— z are all in the interval , where K n is the smallest eigenvalue of P. 

Proof. The proof globally proceeds along the lines of the proof of Lemma 8 in |Sze04) . For x £ X let e x 
denote the characteristic vector of x. Let w\, . . . , w n be the eigenvectors of P with associated eigenvalues 
1 = K\ > . . . > K n > 0. Let v be an arbitrary eigenvector of P z with eigenvalue A. Since P is ergodic, 
\\P Z \\ < 1) therefore A < 1. We show that A > K n . This is obviously true if A = for some k; let us suppose 
that this is not the case. 

Let w be the vector obtained from v by augmenting it with a in the z-coordinate. We express both w 
and e z in the eigenbasis of P: let 10 = X^fc=i IkWk and e z = 53fc=i &kWk- Then wP — Xw + ve z for some real 
number v. Moreover, v ^ 0; otherwise w would have been an eigenvector of P, meaning that A = k^, which 
contradicts our supposition. For k = 1, . . . , k, we have Kklk = A7& + vSk- Since w and e z are orthogonal, we 
also have J2k=i Ik&k — 0. Therefore J2k=i k — a = ^- The statement then follows since the left hand side of 
the above equation would be positive if A were less than n n . □ 

Definition 2.1. For z £ X, the z-hitting time of P, denoted by HT(P, z), is the expected number of steps 
the chain P takes to reach the state z when started in the initial distribution n. 

It is well known that the z-hitting time of P is given by the formula HT(P, z) = i^_ z {I — P_ 2 ) _1 u_ z , 
where u is the all-ones vector. Simple algebra shows that 

it z {i - p- z )- l u. z = y/w^{i-S- z )- x y/wi;, 

where ^Ji^~ z is the entry-wise square root of 7r_ 2 and S- z is the "symmetrized form" SL Z = yTL~7P_ 2 yTTTTT 
of P-z with = diag(jTx)x^z- The matrices P_ z and S- z have the same spectrum since they are similar. 
Let {vj : j < n — 1} be the set of normalized eigenvectors of SL 2 where the eigenvalue of Vj is Xj = cos8j 
with < 9j < 7r/2. By reordering the eigenvalues we can suppose that 1 > Ai > . . . > A n _i > 0. If 
^"K-z = J2j v j v j i s the decomposition of ^tt- z in the eigenbasis of SL Z then the z-hitting time satisfies: 



HT(P,z) = J2j 



v 1 



- X 



j 



When < 9 < it/2 then 1 — 9 2 /2 < cos# < 1 — 9 2 /4. Therefore we can approximate the hitting time with 
another expectation that is very closely related to the analogous quantum notion. More precisely, let H z be 
the random variable which takes the value 1/9 2 with probability z/?, and with probability 1 — J^j v "]- We 
denote the expectation of H z by E[H Z ], Then we have 2E[H Z ] < HT(P, z) < 4E[H Z ]. 
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In the definition of the hitting time the Markov chain is used in a Las Vegas algorithm: we count the 
(expected) number of steps to reach the marked element without error. We can also use the chain as an 
algorithm that reaches the marked element with some probability smaller than 1, leading to a Monte Carlo 
type definition. Technically, to be able to underline the analogies between the classical and quantum notions, 
we define the hitting time with error via the random variable H z . 

Definition 2.2. For z 6 X and for < £ < 1, the e-error z-hitting time of P, denoted by HT e (P, z) is 

defined as 

HT e (P,z) = min {y : Pr[H z > y] < e}. 

Observe that for all z, if e < e' then HT £ /(P, z) < HT £ (P, z). We first show that the use of H z for the 
definition of the Monte Carlo hitting time is indeed justified (proof in Appendix |A")) . For this, let us denote 
by h e (P, z) the smallest integer k such that the probability that the chain does not reach z in the first k 
steps is at most e. 

Theorem 2.2. For all z and e, we have 

h £ (P,z) < 
HT £ (P,z) < 

How much smaller than the Las Vegas hitting time can the Monte Carlo hitting time be? The following 
results state that for some z they are of the same order of magnitude. 

Theorem 2.3. We have the following inequalities between the two notions of hitting time: 

• For all z and e, HT £ (P, z) < ± HT(P, z). 

• There exists z such that for all e < 1/2, HT(P, z) < 4 HT e (P, z). 

Proof. The first statement simply follows from the Markov inequality and from the relation E[P Z ] < 
HT(P, z)/2. For the second statement, let z be an element such that v\ > 1/2. The existence of such 
an element is assured by Lemma 8 in |Sze04j . Then HT(P, z) < TV 4^ 2 /#2 <A/e\< 4HT e (P, z). □ 

3 Quantum hitting times 
3.1 Two notions of quantum walk 

Let U = U2U1 be an abstract search algorithm as in |AKR05j . where [A = I— 2|/i)(/i| for a "target vector" 
with real entries, and U2 is a real unitary matrix with a unique 1-eigenvector |</>o)- Without loss of generality 
we always assume that \<t>o) has real entries. The state \n) is the quantum analogue of the state z which we 
seek in the classical walk P, U2 the analogue of P, and |<^o) the analogue of the stationary distribution it. 

The abstract search algorithm usually starts with state \4>o), and iterates U several times in order to get a 
large deviation from \4>o)- In this paper, we prefer to start with a slightly different initial state. The general 
behavior of the abstract search algorithm remains unchanged by this. We replace the starting state \4>o) by 
|0o) = \(/>o) — ((/>□ I A*) l/J.) , the (unnormalized) projection of the 1-eigenvector |^ ) of U2 on the space orthogonal 
to |/it). This substitution was first considered in |Sze04| . and corresponds to first making a measurement 
according If the measurement outputs we are done. Otherwise we run the abstract search 

algorithm. 

This choice of the initial state is motivated by the results in Section 13.41 which relate quantum hitting 
time to classical hitting time. All other results in this paper remain valid if we keep \4>o). 

Ambainis et al. characterized the spectrum of U in term of the decomposition of \fi) in the eigenvector 
basis of 1/2- One of their results is: 



4 In- ) HT e/2 (P, z), and 



^h e/3 (P,z) 
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Theorem 3.1 ([AKR05]). Let U2 be a unitary matrix with real entries and a unique 1-eigenvector \cj>o}- Let 
\fi) be a unit vector with real entries, and let U\ — I — 2|/i)(/i|. Let U = U$J\. 

• If(<po\(i) = 0, then \<j> ) = \<j) ) and U\(j>o) = \4> Q ) . 

• If (0o I m) 7^ 0; then U has no 1-eigenspace. 

Thus one can use U in order to detect if (0o|a*) 7^ 0. Indeed, in that case, after a certain number T 
of iterations of U on \4>o), the state moves far from the intial state \4>o)- Such a deviation caused by some 
operator V (in our case V = U T , i.e., U iterated T times) is usually detected by using the well known control 
test which requires the use the controlled operator c-V. Namely observe that (H €5 1)(c-V)(H ® 7)|0)|^>) = 
i 1 0} ( I V 7 ) + + — ^IV'))- Therefore a measurement of the first register gives outcome 1 with 

probability HIV') - F|?A)|| 2 /4. 

Szegedy |Sze04j designed a generic method for constructing an abstract search algorithm given a (classical) 
Markov chain. Let P — (p X y) be an ergodic Markov chain over state space X = {1, . . . , n} with stationary 
distribution The time- reversal P* of this chain is defined by equations ^ v p yx — i^xPxy- The chain P is 
reversible if P = P* . 

The quantum analogue of P may be thought of as a walk on the edges of the original Markov chain, 
rather than on its vertices. Thus, its state space is a vector subspace of Ji = C XxX . For a state \ip) £ H, 
let II^, = denote the orthogonal projector onto Span(|^)), and let ref(^) = 2IL/, — Id denote the 

reflection through the line generated by where Id is the identity operator onH. If K. is a subspace of 
H spanned by a set of mutually orthogonal states {\ipi} : i £ I}, then let IL; = ^2 ieI IL^ be the orthogonal 
projector onto IC, and let ref(/C) = 211^ — Id be the reflection through K.. Let A = Spandx)^) : x € -X") 
and B = Span(|p*)|y) : y £ X) be vector subspaces of H, where 

\Px) = J2Vp^v\v) and \P*y) = VPy^\ x )- 

yeX i6Jf 

Define similarly for any z £ X the subspaces A- z = Span(|a;)|p x ) : x £ X \ {z}) and B- z = Span(|p*) \y) : 
A {•-!:•• 

Definition 3.1 ([Szc04 ( MNRS07 ). Let P be an ergodic Markov chain. The unitary operation W(P) = 
rcf(S) • rcf(_4) defined on Ti. is called the quantum analogue of P; and the unitary operation W(P,z) = 
rcf(£>_ z ) • ref(„4_ z ) defined on H is called the quantum analogue o/P_ 2 . 

The unitary operation SWAP is defined by SWAP|x) \y) = \y)\x). When P is reversible, the connection 
between the quantum walk of Szegedy and the quantum walk of Ambainis et al. is made explicit by the 
following fact. 

Fact 3.2. Let z £ X and \fi) = \z)\p z ). Let U 2 = SWAP • ref(^) and Ui = I - 2\fi)(fi\. If P is reversible 
then (U2U1) 2 — W(P, z). In particular, the unitary operators U = U2U1 and W(P,z) are diagonal in the 
same orthonormal basis. 

3.2 Phase estimation and quantum hitting time 

Let U be a unitary matrix with real entries. The potential eigenvalues of U are then 1, —1, and pairs of 
conjugate complex numbers (e IQj , e~ laj ) with < aj < ir, for 1 < j < J, for some J. 

Let be a vector with real entries and of norm at most one. Then uniquely decomposes as 

W = S \w Q )+ 6 3 (\wp + \wJ))+S^\w^), (1) 
i<j<j 

where 8o,S-i,6j are reals, \wq) is a unit eigenvector of U with eigenvalue 1, \w-\) is a unit eigenvector 
with eigenvalue —1, and \wj ), \w~) are unit eigenvectors with respective eigenvalues e laj and e^ IQj , and 
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We now describe a procedure whose purpose is to detect the state has a component orthogonal to 
the 1-eigenspace of U. In the context of the abstract search algorithm, this is equivalent to {4>o\n) ^ 0. The 
idea, similar to the approach of Magniez et al. MNRS07J, is to apply the phase estimation algorithm of 
Kitaev [Kit95j and Cleve et al. [CEMM98] to U. ' 

Theorem 3.3 ( Kit95, CEMM98J). Given an eigenvector \v) of a unitary operator U with eigenvalue e IQ , 
the corresponding phase a £ (— n, n] can be determined with precision A and error probability at most 1/3 
by a circuit Estimate. // \v) is a 1-eigenvector of U , then Estimate determines a = with probability 
1. Moreover, Estimate makes 0(1/A) calls to the controlled operator c-U and its inverse, and it contains 
0((logl/A) 2 ) additional gates. 

Based on the circuit Estimate, we can detect the presence of components orthogonal to the 1-eigenspace 
in an arbitrary state |^>). 



Detect (U, A, e) Input: \ip) 

1. Apply 8(log(l/e)) times the phase estimation circuit 
Estimate for U with precision A to the same state 

2. If at least one of the estimated phases is nonzero, ACCEPT. 
Otherwise REJECT. 



Let QH be the random variable which takes the value 1/a.j with probability 2<5|, the value 1/tt with 
probability S 2 _ 1 , and the value otherwise. Observe that in the following lemma, and in the analysis of all 
our algorithms, the probabilities in fact sum to ||'0I| 2 ; since \ip) is not necessarily normalized, and has norm 
at most 1. 

Lemma 3.4. Assume that Pr[QH > 1/A] < e. Then the procedure Detect (U, A, e) accepts with 
probability ||t/>|| 2 — Sq — 0(e), and moreover with probability if \Sq\ — \\ip\\. In addition, the number of 
applications of c-U is 0(log(l/e)/A). 

Proof. Let us first assume that Estimate can compute the eigenphase of any eigenvector with certainty. 
This assumption is in fact valid when |<5o| = \\ip\\- Then the procedure Detect rejects exactly with probability 

si 

Assume now that Pr[<5-ff > 1/A] < e, and |<$o| < \\ip\\- First observe that Estimate with precision A 
uses 1/A applications of c-U. Then the precision A in Estimate ensures a nonzero approximation of an 
eigenphase ±ctj with probability at least 2/3 provided that Uj > A. By hypothesis, the contribution of these 
eigenphases has squared Euclidean norm 2^^- 5|. The success probability is then amplified to 1 — 0(e) by 
checking that all the 0(log(l/e)) outcomes of Estimate are nonzero. For the special case of eigenphase 0, 
whose contribution has squared Euclidean norm Sq, Estimate gives approximation with probability 1. 

The contribution of the other eigenphases has squared Euclidean norm less than e in the vector 
Therefore the overall acceptance probability is \\tp\\ 2 — Sq — 0(e). □ 

In the case of quantum walk, the above theorem justifies the following definitions of quantum hitting 
times. Let U be some abstract search U2U1, where Ui = I — 2|/i)(/x|, starting from state |</>o) = \(f>o) — ooIa 1 ), 
where ao = (/i|0o). We now set \ip) = \4>q). Again, QH is the random variable which takes the value 1/ctj 
with probability 2<5|, the value 1/ir with probability S_i, and otherwise. 

Definition 3.2. The quantum |/i)-hitting time of U2 is the expectation of QH, that is 

s 2 i 

QHT([/ 2 ») = 2^- + —- 

* — ' OLA 7T 

3 J 

For < e < 1, the quantum e-error \fi) -hitting time of Ui is defined as 

QHT £ ([/ 2 , | M » = minjy : Pv[QH > y] < e}. 



G 



Using Theorem 13. 11 Lemma 13.41 and our definition of quantum hitting time, we directly get: 

Theorem 3.5. For every T > max {1, QHT e (L/2, ^ e procedure Detect {U, 1/T, e) accepts \4>q) with 

probability ||0o|| 2 — 0(e) if ((^oIa*) 7^ 0, and accepts with probability otherwise. Moreover the number of 
applications of c-U is 0(log(l/e) x T). 

If one would like to deal only with normalized states, and to come back to the original starting state \4>o), 
we can encapsulate the projection to the space orthogonal to |/Lt) into our algorithm such as in the following 
main procedure, and deduce its behavior from the above theorem. 



MainDetect([/ 2 , |^), A,e) Input: \ip) 

1. Make a measurement according (\fj,), \/J.) ± ). 

2. If the measurement outputs \fi), ACCEPT. 
Otherwise apply Detect(£7, A, e). 



Corollary 3.6. For every T > max{l, QHT e (£/2, \fj))}, the procedure MainDetect([/2, |^)> ^/T,e) accepts 
\(/)q) with probability 1 — 0(e) if {4>$\p) 7^ 0; ond accepts with probability otherwise. 

When the abstract search is built on the quantum analogue of a reversible Markov chain P and = 
|z)|p z ) for some z, we use the following terminology: 

• The quantum z-hitting time of P is QHT(P,z) = QHT(SWAP • ref(A), \z)\p z )); 

• For < e < 1, the quantum e-error z-hitting time of P is QHT £ (P, z) — QHT £ (SWAP • ref(„4), |z)|p z )). 

With different, more technical arguments, Szegedy proved results similar to Theorem 13.51 albeit with the 
parameter QHT(P, z) for symmetric Markov chains: 

Theorem 3.7 (|Sze04 ). When t is chosen uniformly at random in {1,2,..., [QHT(P, z)~\ }, then the expec- 
tation of the deviation \\(W(P, z))*|<^o) — \<M\\ is f2(||<^ ||)- 

Under certain assumptions, Ambainis et al. [AKR05] have a similar result in terms of the smallest 
eigenphase of U2U1. 

Suppose we wish to not only detect if (fi\4>o) 0, but also to map to Then we are led to a 
procedure different from Detect. One possibility is to try to use U in order to move \ip) to an orthogonal 
state that is closer to \fi). 

Definition 3.3. The {/-rotation of is defined as 6q\wq) + ~ l 1 "^")) + S—i\w-i), where the 

decomposition of \ip) in terms of the orthonormal set of eigenvectors l - "^} °f U is given by Eq. ([I]). 

Fact 3.8. If \if>) is orthogonal to both the 1-eigenspace and the (—l)-eigenspace ofU, then the U -rotation of 
\ip) is orthogonal to \ip). 

This operation can be implemented efficiently by the following procedure with further assumptions on U. 
Namely, we would like U to avoid having any eigenvalue close to — 1. This is naturally the case when we 
consider the quantum analogue of a reversible Markov chain whose eigenvalues are all positive. 



Rotate(C7, A, e) Input: \tp) 

1. Apply 0(log(l/e)) times the phase estimation circuit 
Estimate for U with precision A to the same state \tp)- 

2. If the majority of estimated phases are negative 

Perform a Phase a Flip. 
Otherwise do nothing. 

3. Undo the Phase Estimations of Step[TJ 
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Theorem 3.9. Assume that all eigenvalues e la of U satisfy \a\ < w/2. Then for every T > QHT e ([/ 2 , 

the procedure Rotate([/, 1/T, e) maps \4>o) to a state at Euclidean distance 0(y / e ) from the U -rotation of 

\4>o) ■ Moreover, the number of applications of c-U is 0(log(l/e) x QHT e (£/ 2 , It 1 )))- 

Proof. The proof follows the same argument as in Theorem 13.51 □ 
3.3 Comparison between QHT and QHT e 

We assume henceforth that (4>o\n) 0; otherwise, the problem is trivial — by our definition QHT e (J7 2 , It 1 )) = 
QHT(C/ 2 , = 0. 

The Markov inequality immediately implies the following relationship between these two measures: 
Fact 3.10. For all Ui, It 1 ), ande, 

QHT £ ([/ 2 ») < i QHT({/ 2 »). 

The other direction requires a closer look at the spectral decomposition of U2U1. In this section, we again 
follow the framework of the abstract search algorithm U = U2U1. The eigenvalues of U2 different from 1 are 
either —1 or they appear as complex conjugates e ±i6j , where 9j £ (0, n). For convenience, we assume that 
Q-i = tt, = 6*o < 0i < 62 < . . ., and we always use index j for positive integers. Recall that both and 
the 1-eigenvector \(f>o) of U2 have real entries. Writing \fi) in the eigenspace decomposition of U2 we get 

Im) = aol^ + X^^O^ + l^ 7 )) +o-il^-i>. 
j 

where \4>-i) is a (— l)-eigenvector of t/ 2 , and \4>f) are e ±lSj -eigenvectors of U2 such that a_i and a,j are real, 

\4>-i) has real entries, and \<p~) = \<f>j). Since \<po) has real entries, ao is also real. 

Ambainis et al. [AKR05] (see also Tulsi [Tul08 ) show the following relation between the spectrum of J7 2 
and that of U. 

Lemma 3.11 ( AKR05, Tul08]). The eigenvalues e ±la of the operator U are solutions of the equation: 

«o cot - + J2 S [ cot {-^) + cot \~^) ) ~ ail tan 2 = 

The corresponding unnormalized eigenvectors \w a ) = \fj) + i\w' a ) satisfy (fj,\w' a ) = and 

K) = aoCot||0 o )+^a^cot^^^^+)+cot^^i^j |0T)^ - a_ x tan 1 1</»). 

As in the classical case, we are only able to upper bound QHT by QHT e for some particular target 
states |/it). Therefore, we consider in the following lemma (proof in Appendix |B|) an arbitrary set of or- 
thonormal vectors M = {|/i z )} whose span contains \4>o). In the case of the quantum analogue of a Markov 
chain P as in Definition 13.11 a natural choice for \fi s ) is |z)|pz) for some z in the state space of the Markov 
chain. Recall that \4>q) — \<po) — ao|/i). 

Lemma 3.12. Let M — {\(i z )} be a set of orthonormal vectors with real coefficients in the standard basis, 
such that |0o ) G Span(Af). For every z, let a z be the smallest positive real number a such that e ±1Q are 
eigenvalues of the operator U — U%(I — 2|/i z )(/x z |). Then there exists z such that the length of the projection 
°f |0o) onto the subspace generated by \w az ) and is at least l/v2- 

Corollary 3.13. Let M — {|^z)} be a set of real orthonormal vectors such that |0o) G Span(Af). For all U2 

there exists z such that for all e < 1/2, 

QHT e (£/ 2 > z )) = — . 
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Theorem 3.14. Let M = {\/j, z }} be a set of real orthonormal vectors such that \<f>o) G Span(M). For all U2 
there exists z such that for all e < 1/2, 

QHT([/ 2 > 2 )) < QHT £ ([/ 2 ,| Mz )). 

Proof. This is a consequence of Corollary 13. 131 since QHT([/ 2 , \/J>z)) is by definition at most — . □ 

3.4 Quadratic detection speedup for reversible chains 

Let P be an ergodic Markov chain over state space X = {1, . . . , n}. We further suppose that P is a reversible 
Markov chain with positive eigenvalues, otherwise we simply replace P with 7P + (1 — 7)/, for any 7 < 1/2. 
Let z G X. 

Theorem 3.15. Assume that the eigenvalues of P are all positive. Then we have the following relations: 

• For all z, QHT(P, z) < ^HT(P,z)/2. 

• For all z and e, QHT e (P, z) = ^HT e (P,z). 

Proof. We follow the notation introduced in Sections [2l and [3~T1 Then \ = ^ y/?fx \x)\p x ) , |/x) = |z)|p z ), 
\4>o) — J2xex\{z} \ x ) \Px)- Let ^/tt-z — J2j h 'j v j be the decomposition of ^/7r_ z in the normalized 
eigenbasis of P_ 2 where the eigenvalue of Vj is cosOj, with < 6\ < . . . < Q n -x < tt/2. Let Vj[x] denote 
the ^-coordinate of the vector Vj. We set |^) = J2 x ^z v A x ]\ x )\Px) an d \Q) = J2 y ^z v jiy]\Py)\y)- Then 
\4>o) — Ylj u j\Cj}- F° r every j, the vectors |^-) and \Q) generate an eigenspace of W(P, z) that is also 
generated by two normalized eigenvectors with eigenvalues respectively e 2l8j and e~ 2l0j . This argument is still 
true for SWAP • ref(^4_ z ) when we divide the phases by 2, leading to eigenvalues e j and e -1 ^ (cf. Fact 13.21) . 

~ n— 1 v 2 

Since the length of the projection of |(/f>o) to this eigenspace is z/J, we have QHT(P, z) = J2]=i J~ = E[v / -ff 2 J- 
By the Jensen inequality we get 

QHT(P,z) < y/E[H z ] < y/HT(P,z)/2. 

The second relation above immediately follows from QH 2 — H z . □ 

The same quadratic speed-up as above holds when there are multiple marked elements in the state 
space X. The search algorithm and its analysis are similar and are omitted from this article. 



3.5 On the quadratic speedup threshold 

In this section we consider a broad class of quantum walks defined on undirected graphs. We are able to 
show that for a special case of walks on graphs, the quadratic speedup is tight. 

Let X = {1,2, ...,n}. Our notion of quantum walk can be seen as a walk on the edges of a given 
undirected graph G(X,E). Let H(G) = Span(|xy) : (x,y) £ E) be the Hilbert space that a quantum walk 
on G should preserve. In the rest of this section, we only consider operators and states in H(G) for some 
given G. Observe that SWAP preserves H(G) since G is undirected. 

We introduce a notion of reversibility that is justified by Lemma 13 . 1 71 below. 

Definition 3.4. A quantum walk on an undirected graph G = (X, E) is a unitary C/ 2 = SWAP ■ F defined 
on a subspace ofTL{G), where F is matrix with real entries of the form F = ^2 x ^x\ x )^\ x \ ® ^ x , and where 
U2 has a single 1-eigenvector \4>o)- The quantum walk is reversible when SWAP(|0o)) = \4>o}- 

Observe that the definition implies that | </>o) can be chosen with real entries. This definition of quantum 
walk appears, for example, in the survey paper of Ambainis |Amb03j , see also |San08j . Szegedy considered 
for F x a specific kind of reflection based on Markov chain transition probabilities (see Section 13. ip . 
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Definition 3.5. Let U2 = SWAP ■ F be a quantum walk with unit 1-eigenvector \<f>o) = £V y/W^\x)\<j) x 
where tt x > and \<$> x ) is a unit vector with real entries. Then the classical analogue P = (p xy ) of U2 is 

vy 



defined as p xy = (y\<f) x ^ 2 



Since |<fo) is a 1-eigenvector of SWAP • F we directly state: 
Fact 3.16. Let \i}j x ) = F x \4> x ). Then \4> ) = J2 X V^W)\x) ■ 

Lemma 3.17. The classical analogue P of a quantum walk U2 on G is a Markov chain on G with stationary 
probability distribution it. Moreover, P is reversible if and only if U2 is a reversible quantum walk. 

Proof. First we show that P is a Markov chain on G. For every x, we have 



y y 



^ = V(2/|r) 2 = II^H 2 = l- 



Moreover p xy ^ implies (y\<jf) ^ 0, which implies that (x, y) £ E since \<fio) £ H(G). 

Now we verify that it is a stationary probability distribution. First, tt is a probability distribution since 
\4> x ) for all x £ X and \4>o) are unit vectors. That tt is a stationary probability distribution can be seen from 
the following sequence of equalities which hold for every y £ X: 



/ *K X P X y 



= '^^(xy\<po) 2 by definition of P and \<fio) 

x x 

= ^ir y (x\il) y ) 2 by Fact ETTol 

X 

= 7T y \W)f=7T y . 

For reversibility, observe that we similarly have TT x p xy — (xy\(j>o) 2 8,n&TT y p yx = (yx\<fio) 2 — ((xy|SWAP|0o)) 2 - 
P is reversible when these two expressions are equal for every x, y, which happens precisely when the quantum 
walk U2 is reversible. □ 

Finally, we show that the quadratic speedup is tight in the special case of walks for which all of the 
unitaries F x are reflections. We state the result using the notation above. 

Theorem 3.18. Let U2 = SWAP ■ F be a reversible quantum walk such that F x — 2\<f> x )((f) x \ — L, for all 
x £ X . Then for all z and e, 



QHT e (C/ 2 ,|z)|^)) =QHT £ (P,z) = ^/HT e (P,z). 

Proof. Let U\ = I — 2\z){z\ ® |</> z )(</> z |, for some fixed z. Under the hypothesis of the theorem, (U2U1) 2 
is a product of two reflections ref(.A_ z ) and ref(B_ 2 ), where A- z — Sf>zn{\x)\(j) x ) : x £ X \ {z}) and 
ff_ z = Sp an(|^) \y):y£X\ {z}) = SWAPQA). 

From [Sze04| . we know that the spectrum of (U2U1) 2 is completely defined by the discriminant matrix 
D = A*B, where A = E^J^I^I and B = T^J\<f)\uM- We get that D = ((x|^)(^|y))^ 2 ^. 
The reversibility of 17% guarantees that (xy|7r) = (yx\Tr), which implies that ^flT x ~(y\4> x ) — y/W^{x\(/)y) . Since 

\4> y ) has real entries, we have D = vTLP_ z \/lV , where II = diag(7r x ) x ^ 2 and P_ z is the matrix obtained 
from P by deleting the row and column indexed by z. 

Observe that this discriminant is exactly that of the quantum analogue W(P,z). So W(P,z) and 
(U2U1) 2 are equal up to a basis change which maps J2 X y/W£\x)\p x ) to \<f>o), \z)\p z ) to \z)\cj) z ), and therefore 

T,x?z \f^\x)\Px) to \4> ). □ 
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4 Finding with constant probability 

In this section, we extend a technique devised by Tulsi [Tli108 f° r finding a marked state on the 2D grid in 
time that is the square-root of the classical hitting time. We prove that it may be applied to a larger class 
of Markov chains and target states. The technique may be combined with ideas developed in the earlier 
sections to give an algorithm for the quantum analogue of an arbitrary reversible ergodic Markov chain. 

We use the notation of Section [57X1 In our application, there is no (— l)-eigenvector of U%. Therefore the 
marked state has the following decomposition in an eigenvector basis of V^: 

| M ) = ao |0 o >+ °i(l#> + l^"». ( 2 ) 

l<j<j 

where J is some positive integer. Last, we assume in the rest of this section that (<f>o\fj) 0. 

Lemma 4.1. The vectors and \<po) have the following representation in the basis {\w a }} consisting of 



the eigenvectors of U = U2U1 as given by Lemma \3.11\ = ^ Q ,, 1 p \w a ), and \<fio) = J2 a 0| '""'' -' ' 



Proof. Any vector \tp) may be expressed in the orthogonal basis as \ip) = ^2 a j^jr\w a ). The first 

equation now follows from (w Q |/i) = — i(w£J)|/i) = 1. 

By Lemma f3 . 1 H (4>o\w a ) = (0o|(Im) + i| w a)) = a o + aoicot |-. The second equation follows by combining 
the above with \<fio) — \cf>o) — cjo|m)- □ 

Lemma 4.2. The inner product of the target state \fi) and the U -rotation of \<fro) is 2ao X) Q >o l/^f* P • 

Theorem 13 . 151 shows that the quantum hitting time is bounded by the square-root of the classical hitting 
time when U2 is derived from a reversible Markov chain P, i.e., [7 2 = SWAP ■ ref(^4) in the notation of 
Section I3~2"l This allows for the detection of marked elements (or more generally for checking if (fi\<f>o) 7^ 0) 
and also the creation of the [/-rotation of \<fio) in the stated time bound. However, the overlap of the U- 
rotation of \c/)q) with the target may be o(l). Tulsi |Tul08] discovered a technique, described below, to 
boost this overlap to $7(1) in the case of a quantum walk on the 2D grid. 

Let 6 € [0, 7r/2). Let Rg denote the rotation in C 2 by angle 9: 



Re 



cos 9 ~ sin ( 
sin 9 cos 9 



and let \9) = ijj|0>, and \6^) = rI\1). Define Uf = |0)(0| ® Id + |1)(1| ® U lt and U e 2 = {\9)(9\ ® (-Id) + 
\9 ± )(9 ± \ (g) U 2 ). Then Uf = Id - 2|1)(1| (g) meaning that the modified marked state is Then the 

modified abstract search algorithm becomes: 

T(Ui,U 2 ,0) = U e 2 U{ = (|^|®(-Id) + |^>(^|®L/ 2 )(|0)(0|®Id+|l)(l|®C/ 1 ) (3) 

This is precisely the circuit used by Tulsi: his rotation Rg — Rl in our notation. Tulsi |Tul08| proved that 
the principal eigenvalue of the operator above is closely related to that of the unitary operator U2U1. We 
extend his findings in terms that are more readily used in our context. 

The eigenvalues of C/| are the same as those of U2, except for the addition of the new eigenvalue —1. 
The eigenvectors corresponding to eigenvalues e j are now 16*^)1^). Any state of the form \9)\ip) is a —1 
eigenvector of C/|. 

Fact 4.3. The decomposition of \l)\fi) in the eigenvector basis of f/f * s; 

|1)| M ) = oosfll^) [a \<p )+ "i(\<Pi) + \<t>J))) -sme\0)\n), 

V i<j<j J 

where the coefficients ao,aj are precisely those in Eq. (ty). 
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Lemma 4.4. The eigenvalues e , of the operator T(U\, 1/2,9) are solutions to the equation 
Qq cot ■ 



The corresponding unnormalized eigenvectors \w a ,6) = |1)Ia*) + 'A w 'a e) satisfy (1) lA w 'a e) — and 

l< e ) = cose |^> ^„cot i^o) + ( cot (^-y^) l4) + cot (^y^) 

+ sin0|0> (tea(y) lM) 

Proof. We apply Lemma [3.111 from Section [3~3l with Oq = aocos0, = aj cos 8, a e _ 1 = sin#. Note that U2 
does not have any (— l)-eigenvectors (by assumption), but C/f does. □ 

The target vector in the modified algorithm is |l)|/u)- The start state is chosen to be \4>o,e) = \9 ± )\(f>Q). 
The following are analogous to Lemmata 14. II and 14.21 

Corollary 4.5. The vectors and \4>o.e) have the following representation in the basis {\w a .e}} consisting 

of the eigenvectors of T(f7i, U2, 6) as given by Corollary \4-4\ 

1 ~ cot( — ) 

l 1 )!/^) = / ] |i r^\w a ,e), and \(j) ,e) = (a i cos9) ^\w a ,e)- 

a 9 W Wa M\ a s \\ w ot,8\\ 

Corollary 4.6. The inner product of the target state and the T(Ui, 1/2,0) -rotation of \4>o,e) is given 

by the expression (2a cos6>) J2 a °>o lwlT\\ 2 ' 

We choose for the rest of this section 9 e [0, 7r/2] such that tan 9 = ao cot(ai/2)/10. Let a\ be the 
smallest positive eigenphase of the modified search algorithm T({7i, U2, 9). 

Lemma 14.71 (proof in Appendix |C|) proves that a\ is of the same order as the principal eigenphase a\ of 
the original algorithm U2U1. Lemma 14.81 fproof in Appendix [TJJ is the final piece in our argument. It relates 
the norm of the principal eigenvectors of the modified walk to the norm of the original ones. Both lemmas 
extend corresponding results by Tulsi in the case of the 2D grid. 

Lemma 4.7. There is a unique eigenvalue phase ct\ of the operator T(U%, 1/2,0) in (0, a{\. Moreover, 
cot(a?/2) < 1.01 x cot(ai/2). Therefore ifO < a x < tt/4, then 0.78 xai<flf< a\. 

Lemma 4.8. ||w± ai ,e|| < (3cos#) x ||u>± ai ||. 

We have all the ingredients for the main result of this section. 

Theorem 4.9. Let e > be any constant. Suppose that the squared length of the projection of the state \<f>o) 
onto the principal eigenspace of U2U1 is bounded below by 1 — e. Then, for every T > QHT e (L r 2, |0o))/O-78, 
the procedure Rotate(T(Z7i, U2, &), 1/T, 1/4) maps \4>o,e) to a state with constant overlap with the target 
state \l)\fj). 

Proof. First we prove that T = QHT e (Uf, |l)|/x)) is of the order of QHT £ (C/ 2 , Let I = 2a§(cot 2 ^)/|| w ai || 2 
We know that I > 1 — e. Using Lemma [4.11 we get that QHT e (C/2, — l/ a i- Moreover, by definition, 
T < 1/af. By LemmaO T < 1/(0. 78«i) = QHT £ ([/ 2 , |^>)/0.78. We now get our conclusion by applying 
Corollary Lemmata g7] and [HI and Theorem \%M □ 

We combine the above theorem with Lemma [3. 121 to get our final result. 
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Corollary 4.10. Let P be a state-transitive reversible ergodic Markov chain, and let z be any state. Set 
— \z)\p z )> U\ — I — 2|/i)(/i|, and let Ui be the quantum analogue of P. Then for every e < 1/2 and 
T> QHT £ (J7 2 , \4>o))/0.78, the procedure Rotate(T(t/i, U2, 0), 1/T, 1/4) maps |</>o,e) to a state with constant 
overlap with the target state \l)\fi). 

Proof. The proof is direct once we realize that the conclusions of Lemma 13.121 for a particular element z 
remain valid for any element because of the state-transitivity of P. □ 
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A Proof of Theorem EH 

Proof. For the first statement, we set k = (4 In-) HT E / 2 (P, z), and we denote by Sk the probability that 
z is not reached in the first fc steps. The claim follows if we show that s& < e. It is not hard to see that 
s k = n - z P-z u -z = Ylj (cos 8 j) k . We bound s k from above as follows: 

Sk < E ^ + E ^(cos^) fc 

j:l/0f>HT t/2 (i» j:l/e^<HT e/2 (P,z) 

£ e ' 2+ O-JHT^))'- 

The first summation is at most e/2 by the definition of HT e / 2 (P z), and the bound for the second summation 
is justified since 1/6? < HT e /- 2 (P, z) implies cosOj < 1 — l/(4HT e / 2 (P, z ))- Thus the second term is also at 
most e/2 by the choice of k. 

For the second statement, set k — ^h e / 3 (P, z). Then using the definition of h e /^(P,z) and bounding 
(1 - l/2k) 2k from below by 1/3, we get 

| > •» = E^M) 2fe > E ^{ l -YkT > \ E * 

3 j-.COS 9j>l — J)j " ' j:cos 0j>l-^7 

Now observe that if 1/#| > k, then cos 8j > 1 — Therefore Pi[H z > k] < e, and the statement follows. □ 

B Proof of Lemma 13.121 

Proof. We know from Lemma 13.111 that if 

IM*) = aoM +E a ^ G#) + l^ 7 )) + a -±Mi 

3 

then 

2 a of f a + ?' \ / a; — 6*7 \ \ o r 

a o^ cot 2 = ~ LsHz [ cot [~2~ J +COt \~^2) J +a -i^ tan " 2 ^ 

where < a z < #i < 6* 2 < . . . < n. Since 

cot + cot 



a + 6j} t ^ [a-9jX\ 2 sin a 

cos a — cos ( 



this is equivalent to 



2 a it— » 2 2 sin a 2 a 

2 J: cos a — cos 0,- 2 
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We first claim that there exists z such that |ao,z| 2 < l a M| 2 - Indeed, since |<^o) belongs to Span(M), 

z z z z 

Fix now arbitrarily such a z. To simplify the notation, from now on we refer to a z and to the coefficients 
dj.z without the subscript z. Since all of the terms on the right hand side of Eq. j4|) are positive, it follows 
that 

a 2 sin a ,„. 

cot- > — . (5) 

2 cos a — cos 0i 

Since the right hand side decreases if 9\ is replaced by some 0\ < 8 < 7T, we obtain for every j > 1, 

a 2 sin a 

cot- > -, 6) 

2 cos a — cos Oj 

and 

a 2 sin a a a 

cot— > = 2 tan— > tan—. (7) 

2 cos a — cos 7r 2 2 

We also know that the eigenvectors \w± a ) = + i\w' ±a ) corresponding to the eigenvalues e ±IQ satisfy 



.7 



a cot — \(j) ) + ^ a i ( cot ( 7 -^- ) \(j>f) + cot ( — ) |^ >)-a-i tan 



2 



Let us now define the vector |s) in the two dimensional space generated by |w7± a ) by 

|^ _ \w a ) ~ \W- a ) 

Observe that is orthogonal to \w± a ) and therefore to |s). Then the length of the projection of \(f>o) to 
the subspace is the same as the one of \4>o). This is at least |(s|</>o)|/||s||, which we now bound from below. 
Since the functions tan and cot are odd, we get 

\s) = 2 ao cot||0o)+5>, (cot(^^^ + cot (?^-yj (!</>+> + |^» -2a_i tan | 
and therefore 

,2 



Ml 2 



j 

From (gD, ©, ©, and it follows that 



4aS cot 2 - + 8 > ad 7 — + ^at x tan^ -. 

2 ^ J (cos a - cos 9 j) 2 1 2 



2 

9 o O _ V — ■( 9 Sm Q: 99 

a„ cot — > 2 > a.- — + a , tan 

2 " ^ 3 cos a -cos 6» 2 _1 

and therefore 

IHI 2 < 8a 2 cot 2 |. 
Since (s|</)o) = 2ao cot ^ we can indeed conclude that 

IWo)l > J_ 
INI " V2' 



□ 



15 



C Proof of Lemma 14.7 



Proof. By the definition of ol\ (Lemma 13.111) . 



a cot • 



at -9 



E „ ? cot l^p. = 0. (8) 



Fix any 9 > and define the monotonically decreasing and continuous function / : (0, 9t) i— » R: 

/(jb) = a 2 cot|+^a 2 ( cot (^-^) +COt (^2^)) - tan20tan f- 

We know that lim 5S _ V Q+ /0&) = +°o, hm x ^g 1 _/(x) = — oo, and f(at) < 0. Therefore there is a unique 
a{ £ (0,ai] such that /(a?) = 0. 

From the monotonicity of cot, for x £ (0, ai], we have 

? a i ( cot (^r 1 ) + cot (^r 1 ) ) - E a ^ 2 ( cot + cot ' a ~ 

Using this inequality together with Eq. (jSJ) and the monotonicity of tan, we bound the function / from above 
as follows: 

f(x) > Oq cot — — a 2 , cot — — tan 2 9 tan — 

> Oq cot — — a 2 , cot — — tan 2 tan — 

2 a; 9 a 

> a cot — — 1.01 x cot — , 

where the last inequality comes from the hypothesis < tanf? < ao cot(ai/2)/10. Since /(a?) = 0, we get 
that cot(a?/2) < 1.01 x cot(a x /2). 

We now prove that /(0.78 x ai) < 0, which concludes the proof. In the rest of the proof, we restrict the 
variable x to the interval [0.78 x at, at]- Let (3 be the solution of tan(/3/2) = tan(ai/2)/-\/l-01 in [0, 7r/2). 
Then f(fi)> 0, and therefore a{ > (3. 

Since at > 0, we have tan(ai/2) > ai/2 and 

/3 tan(^) at 

tan — = — ; > 



2 VL01 ~ 2VL01 

Moreover since < (3 < at < 7r/4, we have tan(/3/2) < 2/3/7T, and therefore 

/3 > — ^= x a > 0.78 xai. 
~ 4VL0T 



□ 



D Proof of Lemma 14.81 

Proof. Since the two vectors \w± ai ,e) (respectively, \w± ai )) have the same norm and are orthogonal, it 
suffices to upper bound the following squared norm: 



4cos^<9 



i 2 , cot 2 a?/2 + 2^\ 



sm a. 



cos a\ — cos f 



tan 2 9 tan 2 — 
2 



(9) 
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By Lemma T4. 71 

al cot 2 a\i '2 < 1.01 a 2 , cot 2 a x /2 

and 

tan 2 tan 2 ^- = 0.01 a 2 , cot 2 ^ tan 2 ^- < 0.01 a. 



2 ---o— 2 — 2 - "—o 

by the choice of tan# = ao cot(ai/2)/10 and the monotonicity of cot 2 = 1/ tan 2 on (0,7r/2). Since 

the quantity inside of the square brackets of (J9]) is at most: 



2 a^cot 2 ai/2 + 2^ Qj 2 



/ sinai 



\ cos a\ — cos i 



Hence, 

II WW - w- ai: e\\ 2 < (8cos 2 6») x ||w Q1 - w_ Q1 || 2 . 

□ 
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